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Conversion  Factors,  Non-SI 
to  SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows; 


Multiply 

By 

To  Obtain 

feet 

0.3048 

metres 

inches 

25.4 

millimetres 

miles  (U.S.  statute) 

1.609347 

kilometres 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (mass) 

0.4535924 

kilograms 

tons  (2,000  pounds,  mass) 

907.1847 

kilograms 

1  Introduction 


Background 

Practical  superconducting  magnetic  levitation  (Maglev)  was  pursued  in 
the  United  States  in  the  1960’s.  Development  of  the  concept  continued  for  a 
short  time  in  the  United  States;  but  in  the  1970’s,  Federal  funding  for  Maglev 
vanished  and  development  effectively  ceased.  Foreign  governments,  however, 
continued  development  of  the  Maglev  concept.  Today,  both  Germany  and  Japan 
have  working  prototypes. 

As  a  result  of  the  evolving  foreign  technology  and  increasing 
transportation  needs,  the  United  States’  interest  in  Maglev  was  renewed  in  1990. 
In  December  1990,  the  National  Maglev  Initiative  (NMI)  was  formed  in 
conjunction  with  the  Department  of  Transportation,  the  Corps  of  Engineers,  and 
the  Department  of  Energy.  Its  purpose  was  to  evaluate  the  potential  for  Maglev 
to  improve  intercity  transportation  in  the  United  States  and  to  determine  the 
appropriate  role  for  the  Federal  Government  in  advancing  this  technology.  As 
part  of  its  evaluation,  the  NMI  formed  an  independent  Government  Maglev 
System  Assessment  (GMSA)  team,  composed  of  experts  from  both  the 
Government  and  private  industry.  Its  purpose  was  to  assess  the  technical 
viability  of  potential  U.S.  Maglev  system  concept  designs  (SCD’s).  The 
Structural  Mechanics  Division  (SMD),  Structures  Laboratory  (SL),  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES),  sponsored  through  the  U.S. 
Army  Engineer  Division,  Huntsville  (HND),  was  part  of  this  team  and  provided 
expertise  in  all  areas  related  to  the  guideway  structure  and  dynamic 
vehicle/guideway  interaction  (VGI).  VGl  refers  to  the  dynamic  interaction 
(coupling)  between  two  separate  dynamic  systems:  the  Maglev  vehicle  and  its 
supporting  flexible  guideway. 

As  a  result  of  their  work  as  part  of  the  GMSA  team,  SMD  developed  a 
unique  finite  element  (FE)  analytical  methodology  for  the  smdy  of  VGI  (Ray  et 
al  1995).  The  methodology  has  two  distinct  applications;  it  can  be  used  to 
provide  an  estimate  of  the  vehicle  ride  quality  for  a  given  vehicle  and  guideway 
design  and  to  accurately  predict  the  dynamic  deflections  and  stresses 
experienced  throughout  the  guideway  structure  as  a  result  of  a  vehicle  passage. 
Ride  quality  results  are  necessary  to  design  a  vehicle  suspension  system  and  to 
determine  the  guideway  stiffness  required  to  meet  specific  ride  quality  criteria. 
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Dynamic  structural  analyses  are  necessary  to  produce  safe,  economical,  and 
accurate  guideway  designs. 

As  the  United  States  interest  in  Maglev  grew,  many  companies 
developed  (through  various  funding  sources)  Maglev  SCD’s  of  their  own.  As  a 
means  of  encouraging  this  type  of  development,  HND  provided  technical 
support  to  some  of  these  companies.  One  such  company  was  American  Maglev 
Technology,  Inc.  HND  funded  SMD,  WES,  to  conduct  a  VGI  analysis  of  their 
system  using  their  FE  analytical  methodology.  This  report  presents  a  summary 
of  that  work. 


Objective 

The  objective  of  this  work  was  to  conduct  a  VGI  analysis  of  the 
American  Maglev  System.  Of  particular  interest  in  the  study  was  a  comparison 
of  the  ride  quality  of  the  vehicle  using  two  different  suspension  designs  and  the 
resulting  affects  on  the  guideway  structure. 


Scope 

The  FE  methodology  for  VGI  analyses  as  described  by  Ray  et  al.  (1995) 
was  utilized  for  the  study  herein.  The  details  of  the  American  Maglev  system 
were  obtained  from  Grant  (1994)  and  Davey  (1994  and  1995a,  b).  All  FE 
modeling  was  accomplished  in  two  dimensions  (2-D),  since  sufficient  three- 
dimensional  (3-D)  response  characteristics  (such  as  vehicle  roll  and  yaw 
inertias)  were  not  available  in  the  references. 

In  order  to  study  the  vehicular  response  over  a  significant  length  of 
guideway,  a  series  of  three  two-span  continuous  sections  of  guideway  were 
modeled;  i.e.,  a  total  of  six  spans  at  (90  ft)  each  with  continuous  spans  at  every 
other  support.  Two  different  scenarios  were  studied  for  the  vehicle:  one 
scenario  with  both  an  active  (nonlinear)  suspension  and  a  passive  (linear) 
suspension  acting  together;  and  one  with  the  passive  suspension  by  itself  This 
comparison  was  made  in  order  to  show  the  effectiveness/need  for  the  elaborate 
active  suspension  system.  Separate  analyses  were  made  for  both  of  these 
scenarios  at  vehicle  forward  velocities  of  (100, 200,  and  300  mph).  The  results 
of  these  analyses  are  provided  and  compared  in  the  form  of  data  plots  and  are 
discussed  in  Chapter  4. 
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2  System  Description 


Guideway 

The  design  details  used  in  the  modeling  of  the  guideway  are  described 
by  Grant  (1994)  and  Davey  (1995b).  The  guideway  structure  is  depicted  in 
Figure  1 .  The  superstructure  of  the  guideway  is  a  two-span  continuous 
composite  structure  comprised  of  two  prestressed  concrete  (PSC)  beams  and  one 
cast-in-place  (CIP)  concrete  slab.  The  modulus  of  elasticity  for  the  PSC  beams 
is  4,888  ksi,  while  for  the  CIP  slab  is  equal  to  3,834  ksi.  The  moment  of  inertia 
of  one  PSC  beam  is  equal  to  632,636  in'*.  The  transformed  moment  of  inertia  of 
the  guideway  is  1,999,784  in'*. 

The  “railings”  on  which  the  vehicle  is  supported,  propelled,  and  guided 
are  attached  to  the  top  of  the  guideway  structure.  The  railings  are  a  20.9-in.-high 
by  1.5-in.-wide  sandwich  of  pulltruded  fiberglass  polyvinyl.  They  contain  both 
the  composite  and  active  lift  coils,  as  shown  in  Figure  2. 


Vehicle 

Details  of  the  American  Maglev  SCD  vehicle  were  obtained  from  Davey 
(1994,  1995a,  b).  The  vehicle  is  depicted  in  Figure  3.  It  is  (90  ft)  long  and 
weighs  approximately  40  tons.  As  shown  previously  in  Figure  2,  there  are  four 
“lift  sites”  along  the  length  of  the  vehicle  and  one  lift  site  per  side  at  each  of 
these  locations,  making  a  total  of  eight  lift  sites.  Each  lift  site  consists  of  four 
“C  set”  magnet  pairs  spaced  at  20-in.  centers.  Each  C  set  magnet  pair 
effectively  wraps  around  the  vertical  coils  on  the  guideway,  as  shown  in  Figure 
3.  The  coils  provide  two  mechanisms  for  realizing  lift:  an  active  lift  coil  and  a 
composite  (null  flux)  lift  coil  system.  Together,  they  provide  the  interaction 
between  the  vehicle  and  the  rail. 
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Figure  1 .  Details  of  guideway  structure  for  American  Maglev  SCD 

The  lift  force  associated  with  the  active  lift  coils  is  independent  of  the 
vehicle’s  speed  and  is  defined  by  the  following  polynomial  equation  : 

Fy  (rows)  =  0.01 04y^  -0.1112y^  -  0.0205y  + 1.3998  (1) 

The  variation  of  the  active  lift  coil  force,  Fy,  with  the  vertical  displacement,  y,  is 
shown  in  Figure  4.  The  maximum  lift  will  be  obtained  when  the  active  lift  coils 
are  centered  on  the  shadow  of  the  magnet  faces.  However,  as  seen  in  Figure  4, 
this  corresponds  to  an  unstable  position,  since  the  spring  constant  (the  slope  of 
the  curve  of  Figure  4)  becomes  zero  at  this  displacement  and  assumes  a  negative 
value  with  slight  perturbations.  Therefore,  the  active  lift  coil  has  been  given  an 
initial  offset  of  1 .5  in.,  as  presented  in  Figure  5.  The  factor  that  drives  the  offset 
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is  the  margin  of  safety  desired,  or  the  “overburden.”  As  seen  in  Figure  4,  an 
offset  distance  of  y  =  1 .5  in.  corresponds  to  an  overburden  of  491  lb  (0.245  ton). 
This  reflects  an  1 8-percent  reduction  in  force  from  the  peak  achievable  without 
any  active  current  control  (the  peak  lift  force  equals  2,800  lb  per  C  set  at  35,800 
amps,  and  the  suggested  operating  point  yields  a  lift  force  of  2,308  lb  per  C  set 
at  35,800  amps).  The  spring  constant  is  the  local  slope,  or  the  first 


Figure  2.  Plan  view  of  vehicle’s  lift  sites 
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Figure  5.  Active  lift  coil  offset  for  stable  levitation 

derivative  of  Equation  1,  with  respect  to  the  vertical  displacement,  y.  It  can  be 
expressed  as ; 

k  =  ^  =  0.03 12y^-  0224y  -  0.0205  (2a) 

dy 

At  y  =  1.5  in.,  this  spring  constant  becomes  : 

k  =  -0284  tonslin .  (2b) 

In  a  dynamic  simulation,  the  spring  constant  will  be  equal  to  the  expression 
defined  by  Equation  2a  evaluated  at  the  instantaneous  coil  position,  y.  Figure  6 
shows  a  plot  of  the  stiffness  constant  variation  with  displacement. 
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Figure  6.  Nonlinear  stiffness  variation 

As  mentioned  previously,  there  is  also  a  secondary  mechanism  for 
realizing  lift  known  as  the  composite  (null  flux)  lift  coils.  The  forces  from  these 
coils  are  superimposed  on  those  from  the  active  lift  coils  and  can  be  simulated 
using  a  stationary  winding  excited  at  a  frequency,  co.  This  frequency  is  related 
to  the  velocity  of  the  vehicle  according  to  this  relationship  : 

{in  radfsec)  (3a) 

(in  Hz)  (3  b) 

Upon  substitution  in  Equation  3b,  the  frequencies  corresponding  to  velocities  v 
=  100-,  200-,  and  300  mph,  respectively,  are  /=  43.3-,  86.5-,  and  129.8  Hz, 
respectively.  The  composite  lift  coil  forces  are  therefore  dependent  upon  the 
speed  of  the  vehicle  and  were  calculated  by  Davey  (1995b)  as  : 

F_y  (1 00  mph)  =  589.77y  +  0.808  (4a) 

Fy  (200  mph)  =  644.4y  +  0.7785  (4b) 

F^(300  mph)  =  66A55y  + 1.796  (4c) 

A  plot  of  these  relationships  showing  the  variation  of  the  composite  lift  coil 
force  with  displacement  for  various  speeds  is  presented  in  Figure  7. 


CO  =  A  -v 


/  =  - 
2tc 
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FE  Methodology 

The  FE  methodology  was  utilized  to  simulate  the  dynamic  response  of  the 
Chicle  ^nd  guideway  portions  of  the  SCO  and  the  interaction  between  these  two 
componentsf  known  as  VGI.  The  FE  method  of  structural  analysis  in  general  is 
welUocumented  in  the  literature  and,  thus,  its  theory  is  not  discussed  herem^ 
However  for  purposes  of  modeling  VGI  in  a  Maglev  system,  a  unique  problem 
arises  in  the  application  of  the  FE  method.  The  FE  mesh  for  the  vehicle  system 
(including  its  suspension)  must  be  moved  over  the  FE  mesh  for  Ae  g^^^eway 
structure  to  simulate  the  vehicle’s  forward  motion  across  the  S^^^eway.  This 
was  accomplished  using  interaction  surfaces  referred  to  as 
the  ABAQUS/Standard  code  (Hibbitt,  Karlsson  and  Sorensen,  Inc.  1994).  1  he 
dements  available  in  the  ABAQUS/Standard  elements 
cornerstone  of  the  FE  VGI  methodology.  These  elements  dlow  modeltng  of  the 
interaction  between  deformable  structures  along  sliding  surfaces,  ’^^ere 
separation  and/or  sliding  are  of  finite  amplitude,  and 

surfaces  may  arise.  Slideline  elements  use  a  "master-slave  concept  to  enforce 

the  contact  constmints.  The  surface  on  which  fte  slidelme 
is  the  "slave"  surface.  A  specific  set  of  nodes  is  then  defined  as  the  rnaster 
surface  The  slideline  elements  (or  the  slave  surface)  slide  over  the  slidelme  (o 
master)  surface.  The  nodes  of  the  slideline  elements  are  constrained  not  to 
rSteTnm  the  master  surface.  Generally,  the  master  surface  ,s  chosen  as  the 
See  of  the  stiffer  body  if  the  materials  are  different,  or  as  the  surface  with  the 
ToSer  Sh  The  slave  surface  is  defined  by  attaching  slide  me  mterface 
elements  to  the  surface  of  one  of  the  bodies  in  a  predefined  plane  (that  is 
assumed  to  intersect  the  surface  orthogonally)  and  associating  elements 
with  a  set  of  defined  slideline  nodes  in  the  same  plane.  Relative  "’“tion  alo  g 
theLe  of  interaction  can  be  arbitrarily  large,  but  relative  motions  out  of  the 
plSontlingthe  line  of  interaction  are  neglected,  and  so  must  be  small 
rnmnared  to  typical  element  sizes  on  the  surface. 

:eW:Mincluding  Its  complete  SUSP- 

St)  can  N  molwwith  combinations  of  linear  and  nonlinear  sprtng  and 
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damper  elements.  The  interaction  between  these  two  FE  meshes  can  then  be 
modeled  with  slidelines.  The  slideline  elements  (slave  surfaces)  are  attached  to 
the  vehicle  mesh  at  the  bottom  of  its  suspension  system,  and  the  slideline  nodes 
(master  surfaces)  are  laid  along  the  guideway  mesh  like  the  tracks  for  a  train. 


Guideway  Mesh 

The  FE  mesh  of  the  American  Maglev  guideway  is  depicted  in  Figure  8. 
All  FE  modeling  was  accomplished  in  2-D,  since  sufficient  3-D  response 
characteristics  (such  as  guideway  curvatures,  and  vehicle  roll,  and  yaw  inertial 
capabilities)  were  not  available  on  the  SCD. 


u. 
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Figure  8.  American  Maglev  SCD  guideway  mesh 


In  order  to  study  the  vehicular  response  over  a  significant  length  of  guideway,  a 
series  of  three  two-span  continuous  sections  of  guideway  were  modeled,  i.e.,  a 
total  of  six  spans  at  90  ft  each  with  continuous  spans  at  every  other  support. 
Notice  also  in  Figure  8  that  a  “runup”  and  a  “runout”  span  were  defined  at  each 
end  of  the  guideway  mesh.  These  portions  were  of  the  same  length  as  the 
vehicle  (90  ft)  and  provided  a  starting  point  (prior  to  entering  the  actual  spans 
under  consideration)  and  an  ending  point  for  the  vehicle  in  order  to  get  the 
vehicle  completely  off  the  actual  spans.  The  spans  were  made  very  stiff  not  to 
induce  initial  displacements  into  the  vehicle  prior  to  entering  the  critical  spans. 

The  actual  guideway  spans  were  modeled  using  two-noded  rectangular 
beam  elements,  each  20  in.  long.  They  represented  the  composite  guideway 
section  formed  by  two  PSC  beams  and  a  CIP  slab  (as  discussed  previously). 

They  had  a  width  and  height  of  96-  and  63  in.,  respectively.  These  dimensions 
were  chosen  to  provide  a  transformed  moment  of  inertia  of  1,999,784  in.  ,  with  a 
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corresponding  mass  density  of  9.9  X  10  ’  lb-sec  /in  .  The  mass  of  the 
intermediate  diaphrams  was  neglected.  A  hnear-elastic  material  model  was  used 
to  represent  the  concrete  with  an  elastic  modulus  of  4,888  ksi. 

While  the  actual  guideway  design  is  prestressed,  no  prestressing  or 
reinforcing  was  included  in  the  FE  model  since  these  details  were  not  available, 
and  calculation  of  internal  bending  stresses  was  not  desired  from  these  analyses. 
Conversations  with  the  American  Maglev  designers  indicated  that  die 
prestressing  would  be  used  to  effectively  remove  dead-load  deflections.  No 
mention  was  made  of  using  upward  camber  to  benefit  ride  quaility,  as  has  been 
done  in  the  past  with  other  Maglev  systems  (Lever  1993).  Therefore,  the 
prestressing  effect  was  simulated  in  the  FE  model  through  application  of  an 
upward  uniform  load  equal  to  2.78  k/ft  (equal  to  the  dead  load  of  the 
superstructure)  which  eliminated  dead-load  deflections  and  effectively  provided 
zero  camber  under  dead  load  only. 


Vehicle  Mesh 

The  2-D  FE  mesh  of  the  American  Maglev  vehicle  is  presented  in  Figure 
9.  The  vehicle  body  was  modeled  as  a  series  of  rigidly  connected,  two-noded 
beam  elements  with  dimensions  of  8  x  8  ft  throughout  a  total  length  of  90  ft. 

Since  a  mass  distribution  was  not  specified  in  the  references,  the  total  weight  of 
the  vehicle  (40  tons)  was  distributed  uniformly  over  its  length.  This 
corresponded  to  a  mass  density  of  2.1  x  lO’  lb-sec^  /in.  Vehicle  flexibility  was 
also  undefined.  Since  this  was  considered  to  have  httle  effect  on  the  total 
dynamic  response  of  the  system,  a  very  stiff  modulus  of  elasticity  equal  to 
30,000  ksi  was  used  to  keep  this  mode  of  response  at  a  minimum. 

The  suspension  for  the  vehicle  was  modeled  as  a  series  of  overlapping 
spring  elements  connected  to  the  vehicle  beam  elements  at  their  top  end  and  to 
the  slideline  elements  at  their  bottom  end.  As  seen  in  Figure  9,  the  spnngs  were 
located  along  the  length  of  the  vehicle  at  each  C  set  magnet  pair  location,  as 
previously  shown  in  Figure  2.  Two  overlapping  spring  elements  were  used  to 
represent  each  C  set,  which  experiences  a  combined  lift  force  from  the  composite 
lift  coil  (with  a  linear  force  variation)  and  the  active  lift  coil  (with  a  nonlinear 
force  variation).  A  hnear  spring  was  used  to  represent  the  lift  from  the  composite 
lift  coil,  while  a  nonlinear  spring  was  used  to  represent  that  one  associated  with 
the  active  Uft  coil.  The  springs  were  given  an  unrealistically  long  length  of  4  ft 
to  ensure  a  sufficient  gap  allowance  for  settlement  and  total  dynamic  response. 
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As  shown  in  Figure  7  and  Equations  4a  through  c,  the  force- 
displacement  relationship  for  the  compdsite  lift  coil  is  linear  and  varies  with  the 
vehicle’s  forward  velocity.  The  slope  of  the  lines  in  this  figure,  of  the  first 
deivative  of  Equations  4  a  through  c,  was  used  to  define  the  spring  cons^ts,  k, 
for  the  linear  springs  at  specific  vehicle  velocities.  Since  the  FE  model  is  2-D, 
the  spring  constant  values  were  multiplied  by  2  to  appropriately  represent  the 
right  and  left  side  of  the  C  sets  at  each  location  along  the  length  of  the  vehicle. 


Figure  9.  American  Maglev  SCD  vehicle  mesh 


The  resulting  spring  constants  were  as  follows: 

mph)  =  589.77(2)=  1 17934  Iblin.  (5a) 

ky{20Q  mph)  =  644.40(2)=  1288.80  lb/ in.  (5b) 

*^(300  mph)  =  66435(2)=  1329.10  Ibjm.  (5c) 

The  force-displacement  function  for  the  active  lift  coil  was  previously 
shown  in  Equation  1  and  depicted  graphically  in  Figure  4.  This  function  was 
used  to  define  the  behavior  of  the  nonlinear  springs  of  the  vehicle  mesh.  The 
properties  of  the  nonlinear  springs  were  defined  in  a  numerical  form  inside  the 
ABAQU S/Standard  input  file  (Hibbitt,  Karlsson  and  Sorensen,  Inc.  1994),  which 
is  included  in  Appendix  A.  Careful  review  of  this  segment  of  the  input  file  will 
reveal  that  the  displacement  magnitudes  do  not  correspond  to  those  shown  in 
Figure  4.  The  displacement  values  are  actually  offset  by  3.5  in.,  and  the 
explanation  for  this  is  presented  in  the  following  paragraph. 

As  discussed  previously,  American  Maglev  chose  an  initial  offset 
(overburden)  of  1 .5  in.  on  these  coils  to  prevent  instabilities  in  the  system.  From 
the  curve  shown  in  Figure  4,  it  can  be  observed  that  this  offset  will  correspond  to 
a  lift  force  of  approximately  2,308  lb.  Because  this  is  a  2-D  FE  model,  the 


Chapter  3  Finite  Element  (FE)  Analysis 


13 


nonlinear  spring  force  at  each  C  set  location  will  need  to  be  equal  to  2,308  x  2  = 
4,616  lb  to  properly  simulate  the  total  lift  force.  This  is  the  force  required  by  the 
nonlinear  springs  at  the  beginning  of  each  dynamic  FE  analysis.  The  total 
weight  of  the  vehicle  is  80,000  lb  and,  assuming  the  weight  is  evenly  distributed 
between  all  16  spring  sets,  each  spring  set  (i.e.,  overlapping  linear  and  nonlinear 
springs)  will  support  an  approximate  weight  equal  to  5,000  lb  (calculated  as: 
80,000  lb/ 16  spring  sets  =  5,000  lb).  If  the  nonlinear  springs  support  4,616  lb 
of  the  total  5,000  lb,  then  the  overlapping  linear  springs  must  supplement  the 
supporting  force  by  384  lb  (calculated  as:  5,000  -  4,616  =  384  lb).  The 
deflection,  A ,  of  a  linear  spring  is  defined  as: 


A  =  -  (6a) 

k 

where  F  =  force  in  the  spring,  k  =  spring  constant. 

Therefore,  with  a  required  starting  force  in  the  linear  spring  of  384  lb, 
the  required  initial  displacement  in  the  linear  spring  can  be  calculated  as: 


A  = 


384/Z> 

k 


(6b) 


Using  the  k  values  defined  by  Equations  5a  through  c,  the  required  initial 
displacements  were  thus  equal  to  0.326-,  0.298-,  and  0.289  in.  for  100-,  200-, 
and  300  mph,  respectively.  Thus,  it  can  be  concluded  from  this  discussion  that 
two  specific  initial  conditions  had  to  exist  at  the  beginning  of  the  dynamic  FE 
analysis  to  accurately  represent  the  actual  design  of  the  vehicle:  a  starting  force 
per  nonlinear  spring  equal  to  4,616  lb,  and  an  initial  displacement  of 
approximately  0.30  in.  for  each  linear  spring  to  get  the  required  supplement 
force  of  384  lb.  This  supplemental  force  was  applied  by  artificially  offsetting 
the  displacement  magnitudes  (not  differentials)  of  the  nonlinear  force- 
displacment  relationship  to  the  point  that  a  downward  displacement  of  0.30  in. 
corresponded  to  a  force  of  4,616  lb  in  the  nonlinear  spring. 

In  order  to  verify  that  the  nonlinear  springs  were  properly  defined  in  the 
FE  mesh,  a  simplified  mesh  was  first  analyzed.  It  consisted  of  a  single  lumped 
mass  element  suspended  on  a  single  nonlinear  spring  with  the  same 
characteristics  as  those  used  in  the  actual  vehicle  mesh.  The  lumped  mass  was 
given  a  weight  of  5,000  lb,  which  corresponds  to  1/1 6th  of  the  total  vehicle 
weight,  since  there  are  16  springs  on  the  actual  vehicle  model.  A  dynamic 
analysis  was  performed,  with  a  time  increment  equal  to  0.01  sec  for  a  total  time 
period  of  0.5  sec.  Figure  10  shows  the  resulting  vertical  displacement  history  of 
the  mass  as  it  settled  on  the  spring.  It  can  be  observed  from  this  figure  that  the 
vertical  displacement  followed  a  pattern  similar  to  that  followed  by  the  curve 
shown  in  Figure  4.  This  indicates  that  the  behavior  of  the  nonlinear  spring 
element  was  properly  defined. 
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Siidelines 


As  described  in  a  previous  section,  the  *SLIDELINE  option  in  the 
ABAQUS/Standard  code  (Hibbitt,  Karlsson  and  Sorensen,  Inc.  1994)  was  used 
to  provide  the  sliding  interaction  between  the  vehicle  and  the  guideway  mesh. 
The  slideline  elements  were  attached  between  adjacent  nodes  corresponding  to 
the  bottom  ends  of  the  vehicle  spring  elements.  Their  nodal  connectivity  was  set 
as  a  hinge-type  connection  to  allow  independent  vertical  responses  at  each  of  the 
spring  elements,  while  still  providing  the  sliding  interaction  along  the  guideway 
(i.e.,  the  slave  surface).  A  continuous  slideline  surface  (i.e.,  the  master  surface) 
was  defined  along  the  length  of  the  guideway  mesh  and  between  the  guideway 
nodes. 


Figure  10.  Variation  of  the  active  lift  force  with  displacement 
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Dynamic  Analysis 


The  ABAQUS/Standard  FE  code  (Hibbitt,  Karlsson  and  Sorensen,  Inc. 
1994)  was  used  for  the  dynamic  analysis  of  the  American  Maglev  mesh.  A 
typical  analysis  input  file  is  included  in  Appendix  A.  The  vehicle  was 
effectively  moved  across  the  guideway  mesh  at  a  constant  forward  velocity  by 
applying  a  velocity  boundary  condition  to  all  of  the  nodes  associated  with  the 
vehicle  mesh.  Analyses  were  made  for  vehicle  velocities  of  100-,  200-,  and  300 
mph.  Separate  analyses  were  also  made  using  only  the  linear  spring  suspension 
(i.e.,  the  composite  coil  lift)  for  the  same  velocities.  For  these  analyses,  the 
linear  spring  constants  in  Equations  5  a  through  c  were  multiplied  by  a  factor 
equal  to  12.3 115,  according  to  Davey  (1994). 

Each  FE  analysis  was  run  for  the  total  time  required  to  move  the  vehicle 
completely  from  one  end  of  the  guideway  to  the  other  (i.e.,  over  seven  spans, 
including  the  runout  section).  These  times  corresponded  to  4.36,  2.16,  and  1.43 
sec  for  speeds  of  100,  200,  and  300  mph,  respectively.  Time-steps  of  0.01 15, 
0.0057,  and  0.0038  sec,  respectively,  were  utilized  for  each  of  these  analyses. 
These  time  steps  were  chosen  to  ensure  that  each  guideway  element  received  at 
least  one  force  interaction  from  each  spring  as  the  vehicle  moved  over  it.  This 
was  calculated  with  the  following  equation: 


,  .  .  Guideway  element  length 

Analysis  time  -  step  = - 

Vehicle  velocity 


(7) 
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4  Analytical  Results 


The  results  obtained  from  the  dynamic  analyses  are  presented  in  Figures 
1 1  through  30.  Each  analysis  was  performed  for  two  different  scenarios:  the 
vehicle  when  supported  only  on  its  composite  lift  coil  suspension  system 
(referred  to  herein  as  “Linear  Springs  Only”),  and  with  the  vehicle  supported  by 
both  the  composite  lift  and  active  lift  coil  systems  (referred  to  herein  as  “Both 
Springs”).  All  the  results  presented  herein  were  nondimensionalized  according 
to  the  position  of  the  front  of  the  vehicle  (i.e.,  head  position)  with  respect  to  each 
individual  90-ft  span  length  of  guideway.  This  was  done  by  multiplying  the 
time-axis  (x-axis)  by  the  vehicle’s  velocity  and  dividing  this  value  by  the  span 
length  (90  ft  or  1,080  in.).  Note  that  data  after  the  vehicle  entered  the  runout 
span  are  not  shown,  since  it  was  given  an  unrealistically  large  stiffness.  Also 
note  that  many  of  the  data  plots  stop  prematurely  (i.e.,  prior  to  the  vehicle 
reaching  the  end  of  span  6).  This  situation  occurred  mostly  with  the  linear- 
springs-only  scenario  and  worsened  as  the  vehicle  velocities  increased.  This  was 
most  likely  due  to  too  rapid  a  change  in  displacement  of  the  vehicle  (shown  in 
later  plots)  which  caused  a  convergence  tolerance  problem  within  ABAQUS. 
While  this  was  an  internal  convergence  problem,  it  indicates  possible 
insufficiencies  with  the  suspension  system  that  utilizes  linear  springs  only. 

Figures  1 1  through  16  show  a  comparison  of  the  vertical  displacement 
time-histories  for  the  guideway  midspan  locations  for  spans  1,  2,  and  5.  These 
displacements  were  obtained  for  vehicle  velocities  of  1 00  and  300  mph.  As  can 
be  seen,  the  300  mph  vehicle  velocity  consistently  caused  the  highest 
displacements  of  the  guideway.  As  expected,  the  suspension  system  of  the 
vehicle  made  little  difference  as  to  guideway  response.  For  both  velocities,  the 
peak  responses  consistently  occurred  only  after  the  vehicle  had  completely 
covered  the  span  of  concern.  The  maximum  downward  response  for  the  linear 
springs  only  system  was  0.17  in.  on  span  2  and  for  the  system  with  both  springs, 
was  0.13  in.  on  span  2.  For  spans  of  90  ft  (or  1,080  in.),  these  correspond  to 
very  high  span-deflection  ratios  of  L/6,400  and  L/8,300,  respectively.  For  a 
static  midspan  deflection  of  0.095  in.  (two-span  beam  with  a  uniformly  applied 
vehicle  load  on  one  span),  these  correspond  to  dynamic  load  factors  (i.e.,  ratio  of 
dynamic  to  static  deflection)  for  the  guideway  of  1.79  and  1.37,  respectively. 

All  of  these  values  are  well  within  previously  established  design  criteria  for 
other  U.S.  Maglev  systems  (Lever  1993). 
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The  vertical  displacements  of  the  vehicle  as  it  traversed  the  guideway 
spans  are  presented  in  Figures  17  through  20.  The  displacements  shown  are 
those  of  the  vehicle  nodes  at  the  attachment  points  for  Spring  1  (denoted  as 
“Vehicle’s  Front”)  and  Spring  16  (denoted  as  “Vehicle’s  Rear”).  These 
locations  are  shown  in  Figure  9.  From  these  plots  it  can  be  seen  that  the  300- 
mph  velocity  caused  the  greatest  deflections  for  the  case  of  the  vehicle  on  linear 
springs  only.  However,  when  both  suspensions  were  included,  the  exact 
opposite  was  true;  the  100-mph  velocity  caused  the  highest  response  and  the 
300-mph  velocity  caused  the  lowest.  Note  that  with  both  springs  and  at  300 
mph,  the  vehicle  displacement  continued  to  increase  with  each  successive  span. 
This  result  shows  the  importance  of  multiple-span  analyses  as  conducted  herein 
and  indicates  the  potential  need  for  damping  of  the  vehicle’s  suspension  system. 

The  vertical  accelerations  corresponding  to  the  same  locations  as  the 
above  displacements  are  shown  in  Figures  21  through  24.  For  the  case  of  the 
linear  springs  only,  the  300-mph  velocity  produced  considerably  higher 
accelerations  in  the  0.9-g  range  as  compared  to  values  less  than  0. 1  g  for  the 
100-mph  velocity.  In  addition,  the  increasing  trend  of  the  300-mph  curve 
indicates  that  accelerations  over  successive  spans  could  be  even  greater.  The 
accelerations  for  the  vehicle  with  both  springs  were  much  lower,  with  the  worst 
accelerations  of  approximately  0.12  g  occurring  at  the  200-mph  speed.  The 
frequency  of  all  of  these  accelerations  was  in  the  1-  to  5-Hz  range.  While  these 
data  show  that  the  addition  of  the  active  suspension  made  a  substantial 
improvement  to  the  vehicular  ride  quality,  the  accelerations  were  still  somewhat 
high  compared  to  those  for  other  SCO’s  as  studied  by  Lever  (1993).  They 
showed  that  for  accelerations  in  the  1-  to  5-Hz  range  (as  with  the  data  herein), 
the  ISO  1-hour  reduced-comfort  level  is  between  0.05  and  0.1  g.  The 
accelerations  calculated  herein  were  at  or  above  these  limits. 

Figures  25  and  26  show  the  variation  in  the  nonlinear  spring  forces  for 
the  front  and  the  rear  part  of  the  vehicle  as  it  traveled  through  the  guideway.  As 
desired,  the  force  in  each  spring  started  at  4,616  lb  to  represent  the  1 .5-in.  offset 
in  the  nonlinear  springs.  Throughout  the  guideway  traversal,  the  force  variations 
stayed  within  a  200-lb  range  (or  100-lb  per  C  set)  for  the  300-mph  vehicle  speed. 
The  variations  were  as  much  as  450  lb  (or  225-lb  per  C  set)  for  the  1 00-mph 
speed.  These  variations  were  all  well  within  the  design  force  overburden  or 
safety  margin  (see  previous  discussion  of  the  suspension  system)  on  these 
springs  of  491 -lb  per  C  set. 

Figures  27  through  30  show  the  change  in  length  associated  with  the 
front  and  rear  springs,  which  is  representative  of  the  vertical  gap  variation 
between  the  vehicle  and  the  guideway  throughout  the  vehicle  traversal.  Because 
of  its  null-flux  type  of  suspension  system,  this  is  not  as  important  for  this 
system.  The  horizontal  gap  variation  would  be  more  important,  but  since  this 
was  only  a  2-D  analysis,  these  data  of  course  were  not  obtained. 
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Figure  25.  Spring  force  at  vehicle's  front  (both  springs) 


Figure  26.  Spring  force  at  vehicle’s  rear  (both  springs) 
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Figure  27.  Magnetic  gap  variation  at  vehicle’s  front  (linear  spring  only) 


Figure  28.  Magnetic  gap  variation  at  vehicle's  front  (both  springs) 
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5  Conclusions  and 
Recommendations 


The  FE  methodology  was  successfully  used  to  study  d3mamic  VGI  for 
the  American  Maglev  system.  Modeling  of  the  active  suspension  system  for  the 
vehicle  proved  to  be  especially  challenging,  but  the  resulting  FE  model  proved 
to  represent  the  system  design  quite  well. 

The  analytical  results  showed  the  guideway  superstructure  to  be  very 
stiff  and  to  respond  well  under  the  dynamic  vehicular  loadings.  The  300-mph 
vehicle  speed  consistently  produced  the  greatest  midspan  deflections  on  the 
guideway  and  indicated  that  a  dynamic  load  factor  of  approximately  1 .4  should 
be  used  for  future  guideway  design  applications  with  this  system  (assuming  that 
both  the  active  and  linear  suspensions  are  used  on  the  vehicle).  While  no  stress 
computations  were  made,  the  low  guideway  deflections  indicate  that  high 
stresses  should  not  be  a  problem  in  the  guideway. 

As  would  be  suspected,  the  suspension  system  using  the  combination 
active  and  linear  springs  was  clearly  better  from  the  ride  quality  standpoint. 
Vehicular  accelerations  were  less  by  almost  a  factor  of  10.  However,  while 
lower  than  the  system  with  the  linear  springs  only,  the  accelerations  were  still 
slightly  above  the  ISO  1-hour  reduced-comfort  level  as  discussed  by  Lever 
(1993).  If  desired,  ride  quality  could  possibly  be  improved  through  the 
introduction  of  guideway  precamber,  a  stiffer  guideway,  or  improved  force 
control  characteristics  for  the  active  suspension.  It  should  be  remembered 
however  that  this  was  by  no  means  a  comprehensive  ride  quality  analysis,  and 
the  comments  above  should  only  be  used  as  an  initial  assessment.  A 
comprehensive  ride  quality  analysis  would  also  include  random  guideway 
roughness  and  irregularity  and  would  be  run  over  a  larger  number  of  guideway 
span  lengths. 

While  this  was  by  no  means  an  exhaustive  analysis  of  the  American 
Maglev  system,  it  has  shown  the  system  to  be  viable  from  the  VGI  standpoint. 
Future  studies  in  this  area  should  concentrate  on  studying  3-D  aspects  of  the 
dyamic  response  of  the  system,  with  emphasis  on  the  vehicle.  Since  this  system 
uses  a  null-flux  levitation  concept  where  the  vehicle  magnets  wrap  around  the 
guideway  coils,  horizontal  movements  of  the  vehicle  and  horizontal  magnetic 
gap  variations  will  be  very  important.  A  3-D  analysis  should  consider  wind, 
guidance,  and  centrifugal  force  effects  on  the  vehicle  and  guideway.  The  effect 
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of  these  forces  on  the  fiberglass  railings  (containing  the  guideway  coils)  should 
also  be  carefully  studied. 
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215, 

216, 

217, 

222, 

223, 

224, 

225, 

230, 

231, 

232, 

233, 

238, 

239, 

240, 

241, 

246, 

247, 

248, 

249, 

254, 

255, 

256, 

257, 

262, 

263, 

264, 

265, 

270, 

271, 

272, 

273, 

278, 

279, 

280, 

281, 

286, 

287, 

288, 

289, 

294, 

295, 

296, 

297, 

302, 

303, 

304, 

305, 

310, 

311, 

312, 

313, 

318, 

319, 

320, 

321, 

335, 

336, 

337, 

338 

343, 

344, 

345, 

346 

351, 

352, 

353, 

354 

359, 

360, 

361, 

362 

341, 

345, 

348, 

352 

367, 

371, 

374, 

377 
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*ELSET,  ELSET=PASSUSP 


378, 

379, 

380, 

381, 

382, 

383, 

384, 

385, 

386, 

387, 

388, 

389, 

390, 

391, 

392, 

393, 

★  * 

**  Carbody 

if  "k  _  —  •  . 

— 

*NSET, 

NSET=CARBODY 

330, 

331, 

333, 

335, 

338, 

339, 

341, 

343, 

346, 

347, 

349, 

351, 

354, 

355, 

357, 

359, 

362, 

363, 

★  ★ 

*ELSET, 

ELSET= 

CARBODY 

325, 

327, 

329, 

333, 

337, 

340, 

342, 

346, 

350, 

353, 

355, 

359, 

363, 

366, 

368, 

372, 

376, 

** 

**  Set  the  EQUATION  multi-point-constraints  : 


♦EQUATION 

2, 

331, 

2, 

1, 

1.0, 

332, 

1, 

-1.0 

333, 

2, 

1, 

1.0, 

334, 

1, 

1 

H 

O 

335, 

2, 

1, 

1.0, 

336, 

1, 

O 

H 

1 

338, 

★  ★ 

2, 

1, 

1.0, 

337, 

1, 

-1.0 

339, 

2, 

1, 

1.0, 

340, 

1, 

1 

H 

O 

341, 

2, 

1, 

o 

H 

342, 

1, 

-1.0 

343, 

2, 

1, 

H 

O 

344, 

1, 

O 

H 

1 

346, 

** 

2, 

1, 

O 

H 

345, 

1, 

O 

H 

1 

347, 

2, 

1, 

H 

O 

348, 

1, 

-1.0 

349, 

2, 

1, 

1.0, 

350, 

1, 

O 

H 

1 

351, 

2, 

1, 

1.0, 

352, 

1, 

-1.0 
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1 


-1.0 


354,  1,  1.0,  353, 
** 


2, 

355, 

1. 

1.0, 

356,  1, 

o 

H 

1 

2, 

357, 

1, 

1.0, 

358,  1, 

1 

O 

2, 

359, 

1, 

1.0, 

360,  1, 

O 

H 

2, 

362, 

1, 

1.0, 

361,  1, 

1 

H 

O 

** 


**  Slide  Line  Magnets 
_ _ 

*NSET,  NSET=MAGNETS 


332,  334,  336, 
348,  350,  352, 

337, 

353, 

340, 

356, 

342, 

358, 

344, 

360, 

345 

361 

** 

♦ELSET,  ELSET=MAGNETS 

331,  336,  344, 

349, 

357, 

362, 

370, 

375 

♦  ★ 


**  Slide  Lines 

♦SLIDE  LINE,  ELSET=MAGNETS ,  TYPE=LINEAR,  GENERATE, 
SMOOTH=0.2 


329, 

329, 

0 

328, 

274, 

-1 

273, 

220, 

-1 

219, 

165, 

-1 

164, 

111, 

-1 

110, 

56, 

-1 

55, 

2, 

-1 

1, 

1, 

0 

*  * 

**  Element  Properties 


★  ★ 

**  Guideway 


♦BEAM  SECTION,  SECT=RECT,  MATERIAL=GUIDEMAT ,  ELSET=GUIDEWAY 
96 . , 63  . 

0.,  0.,  -1. 
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★  ★ 


**  End  of  Guideway 


♦BEAM  SECTION,  SECT=RECT,  MATERIAL=ENDMAT,  ELSET=ENDTRACK 
24., 8. 

0.,  0.,  -1. 

**  Suspensions 


**  (ACTIVE  LIFT  COILS)  -->  All  Stiff.  K  =  Ko  *  2 


♦SPRING,  ELSET=ACTSUSP,  NONLINEAR 


/  t 

-5599.2,  -1.8 

-5586.5936,  -1.7 

-5565.3408,  -1.6 

-5535.6912,  -1.5 

-5497.8944,  -1.4 

-5452.2,  -1.3 

-5398.8576,  -1.2 

-5338.1168,  -1.1 

-5270.2272,  -1.0 

-5195.4384,  -0.9 

-5114.0,  -0.8 

-5026.1616,  -0.7 

-4932.1728,  -0.6 

-4832.2832,  -0.5 

-4726.7424,  -0.4 

-4615.8,  -0.3 

-4499.7056,  -0.2 

-4378.7088,  -0.1 

-4253.0592,  0.0 

-4123.0064,  0.1 

-3988.8,  0.2 

-3850.6896,  0.3 

-3708.9248,  0.4 

-3563.7552,  0.5 

-3415.4304,  0.6 

-3264.2,  0.7 

-3110.3136,  0.8 

-2954.0208,  0.9 

-2795.5712,  1.0 

-2635.2144,  1.1 

-2473.2,  1.2 

-2309.7776,  1.3 

-2145.1968,  1.4 
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-1979.7072, 

1.5 

-1813.5584, 

1.6 

-1647.0, 

1.7 

-1480.2816, 

1.8 

-1313.6528, 

1.9 

-1147.3632, 

2.0 

-981.6624, 

2.1 

-816.8, 

2.2 

-653.0256, 

2.3 

-490.5888, 

2.4 

-329.7392, 

2.5 

-170.7264, 

2.6 

-13.8, 

2.7 

140.7904, 

2.8 

292.7952, 

2.9 

441.9648, 

3.0 

588.0496, 

3.1 

730.8, 

3.2 

869.9664, 

3.3 

1005.2992, 

3.4 

1136.5488, 

3.5 

1263.4656, 

3.6 

1385.8, 

3.7 

1503.3024, 

3.8 

1615.7232, 

3 . 9 

1722.8128, 

4.0 

1824.3216, 

4.1 

1920.0, 

4.2 

2302.2, 

4 . 3 

2501.2, 

4.4 

★  * 

**  (COMPOSITE  COILS)  -->  Stiff  Ko=664.55  Ib/in  *  2  sides 
**  _ =============  =>  1329.1  Ib/in 


♦SPRING,  ELSET=PASSUSP 


1329.1, 

★  ★ 

**  Carbody 
★  ★  — 

♦BEAM  SECTION,  SECT=RECT,  MATERIAL=VEHICMAT ,  ELSET=CARBODY 
96 . , 96  . 

0.,  0.,  -1. 

★  ★ 

♦♦  Slide  Line  Magnets 
★  *  —  —  —  —  —  -  —  —  —  —  —  “  — 
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♦INTERFACE,  ELSET=MAGNETS 

1.  ,  0. ,0. ,1. 

♦FRICTION 
0.3 
★  * 

*********************♦♦♦♦*♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦ 
★  ★ 

♦♦  Material  Properties 

************************************************************* 

** 

**  Guideway  Material 

"kif  — 

♦MATERIAL,  NAME=GUIDEMAT 

♦DENSITY 

0.0000992063, 

♦ELASTIC,  TYPE=ISO 
4888. E+03,  0.15 

♦TRANSVERSE  SHEAR  STIFFNESS 
2.17E+06 
★  * 

♦♦  End  Track  Material 

'k  it  — 

♦MATERIAL,  NAME=ENDMAT 

♦DENSITY 

lE-6, 

♦ELASTIC,  TYPE=ISO 
30.E6,  0.3 

** 

♦♦  Carbody 

ie  ie 

♦MATERIAL,  NAME=VEHICMAT 
♦DENSITY 
0 . 000020797 
♦ELASTIC,  TYPE=ISO 
30 .E+10 
*  * 

************************************************************ 
★  ★ 

**  Boundary  Conditions 

★  * 

************************************************************ 
★  * 
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**  Hinges 
** _ 

♦BOUNDARY 
56,  1,2 
165,  1,2 
274,  1,2 

*  * 

**  Rollers 
** 

♦BOUNDARY 

2,  2,2 

★  ★ 

110,  2,2 
111,  2,2 

*  ★ 

219,  2,2 

220,  2,2 

*  * 

328,  2,2 

*  * 

♦♦  Fixed- Ends 
♦* 

♦BOUNDARY 


1, 

1,2 

329, 

1,2 

330, 

1,1 

332, 

1,1 

** 

SPRINGS ,3,3 


★  * 

**  Output  Request 

★  ★ 

★  ★ 

♦*  Applying  Initial  Displacement  on  Springs 
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*STEP,  AMPLITUDE=STEP 
♦STATIC 
0.05,1.0,0.05 
♦BOUNDARY 
CARBODY,2,2, -0.289 
** 

♦♦  Self  Weight  of  the  Vehicle  and  Guideway 
★  *  _ 

♦DLOAD 

CARBODY,  GRAV,  386.21,  0.,-l.,0. 

GUIDEWAY,  GRAV,  386.21,  0.,-l.,0. 

ir  ★ 

♦♦  Distributed  Load  Along  Guideway  (to  consider  pre¬ 
stressing  effect) 

**  _ _ 

♦DLOAD 

GUIDEWAY, PY, 231 . 67 
** 

♦NODE  FILE 

U 

RF 

♦EL  FILE,  ELSET=PASSUSP 
S 

♦EL  FILE,  ELSET=ACTSUSP 
S 

♦NODE  PRINT,  F=0 
♦EL  PRINT,  F=0 
♦END  STEP 
★  ★ 

♦♦  Equilibrating  the  system 
★  ★  —  —  —  —  —  -  —  —  —  —  —  —  —  —  - 

♦STEP,  AMPLITUDE=STEP 
♦STATIC 

♦BOUNDARY,  OP=NEW 


56, 

1,2 

165, 

1,2 

274, 

1,2 

2, 

2,2 

110, 

2,2 

111, 

2,2 

219, 

2,2 

220, 

2,2 

328, 

2,2 

1, 

1,2 

329, 

1,2 
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330,  1,1 
332,  1,1 
SPRINGS,  3,3 
♦NODE  FILE 
U 
RF 

*EL  FILE,  ELSET=PASSUSP 
S 

*EL  FILE,  ELSET=ACTSUSP 
S 

♦NODE  PRINT,  F=0 
♦EL  PRINT,  F=0 
♦END  STEP 
** 

♦♦  Dynamic  Analysis  of  the  System  (Velocity  =  300  mph 
[5275.58  in/sec]) 

**■  - - -  -  -  -  - - - - 

♦STEP,  INC=500 
♦DYNAMIC , DIRECT , NOHAF 
0.003791,1.4330 
♦RESTART,  WRITE,  FREQ=3 
♦PRINT,  CONTACT=yES 
♦NSET,  NSET=MAGLEV 
MAGNETS, 330, 363 
♦BOUNDARY , OP=NEW 
56,  1,2 
165,  1,2 
274,  1,2 
2,  2,2 
110,  2,2 
111,  2,2 

219,  2,2 

220,  2,2 

328,  2,2 

1,  1,2 

329,  1,2 
SPRINGS, 3, 3 

♦BOUNDARY,  TYPE=VELOCITY,  OP=NEW 

MAGLEV,1,1,5275.58 

♦NODE  FILE,  FREQ=3 

U 

V 

A 

RF 

♦EL  FILE,  FREQ=3,  ELSET=PASSUSP 
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s 

*EL  FILE,  FREQ=3,  ELSET=ACTSUSP 
S 

♦FILE  FORMAT,  ASCII 
♦NODE  PRINT,  F=0 
♦EL  PRINT,  F=0 
♦END  STEP 
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